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The predissociation dynamics of lithium iodide (Lil) in the first excited A-state is investigated 
for molecules in the gas phase and embedded in helium nanodroplets, using femtosecond pump- 
probe photoionization spectroscopy. In the gas phase, the transient Li + and Lil + ion signals feature 
damped oscillations due to the excitation and decay of a vibrational wave packet. Based on high-level 
ab initio calculations of the electronic structure of Lil and simulations of the wave packet dynamics, 
the exponential signal decay is found to result from predissociation predominantly at the lowest 
avoided X-A potential curve crossing, for which we infer a coupling constant Vxa = 650(20) cm -1 . 
The lack of a pump-probe delay dependence for the case of Lil embedded in helium nanodroplets 
indicates fast droplet-induced relaxation of the vibrational excitation. 


I. INTRODUCTION 

Three decades ago, the dream of controlling chemi¬ 
cal reactions with laser pulses has initiated the field of 
coherent control [1-3]. The idea is to selectively con¬ 
trol the outcome of photo-induced reactions - formation 
of a chemical bond or creation of photofragments - by 
changing the parameters of the laser pulses. Femtosec¬ 
ond pump-probe spectroscopy has been instrumental in 
bringing this dream closer to reality, terming the notion 
“femtochemistry” [4-7]. 

Alkali halide diatomics were among the first molecules 
to be studied by femtosecond spectroscopy in the gas 
phase. Their lowest electronic excitations mostly lie in 
the visible or ultraviolet spectral range and are there¬ 
fore well accessible by femtosecond lasers. Alkali halides 
are special in that their chemical bond has ionic char¬ 
acter even in the electronic groundstate. This charac¬ 
ter results from the extremely differing electron affini¬ 
ties of the constituent atoms. As a consequence, even 
the lowest excited states are subject to curve crossings 
with the ion-pair state potential curve at intermediate 
distances between the nuclei. These curve crossings give 
rise to non-adiabatic couplings and to predissociation of 
the photo-excited molecules. 

In pioneering time-resolved experiments, Zewail and 
coworkers have visualized the predissociation dynamics 
of sodium iodide (Nal) and other alkali halides as the 
repeated passage of a coherent vibrational wave packet 
across the avoided curve crossing of the groundstate and 
the first excited state [8-10]. At every passage the wave 
packet splits up into one fraction which remains in a 
bound state and continues to vibrate, and another frac¬ 
tion which follows the dissociative branch of the poten- 
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tial, producing unbound Na and I fragments [8, 9, 11, 12], 

Such textbook quantum molecular dynamics can typ¬ 
ically only be observed in gas phase experiments where 
the molecules are isolated from their environment. In 
the condensed phase, the coherent motion induced by a 
femtosecond laser pulse is most often subject to very fast 
decay due to relaxation and dephasing processes. A not- 
icable exception present diatomic alkali metal molecules 
attached to helium nanodroplets [13-17]. Helium (He) 
droplets are quantum fluid clusters which feature extraor¬ 
dinary properties such as microscopic superfluidity [18], 
a critical velocity for the undamped motion of embed¬ 
ded molecules [19], and quantized states of angular mo¬ 
mentum (vortices) [20]. The dynamical response of He 
droplets can be probed by the wave packet motion of 
impurities such as alkali dimers. The vibrational dynam¬ 
ics of these molecules was found to be only weakly per¬ 
turbed by the droplets. This can be explained by a weak 
coupling of alkali metal atoms and small molecules to 
helium droplets due to their position at the droplet sur¬ 
face. In contrast, recent studies of the impulsive align¬ 
ment of molecules embedded in the droplet interior have 
evidenced that the rotational wave packet dynamics is 
significantly slowed down and rotational recurrences are 
completely suppressed [21]. 

Here, we present a femtosecond pump-probe photoion¬ 
ization study of free Lil and of Lil embedded in helium 
nanodroplets. Lithium iodide, as all other molecules 
which do not purely consist of alkali or alkaline earth met¬ 
als, are immersed in the interior of He droplets. Therefore 
strong perturbations of vibrational dynamics by the He 
environment may be expected. We start by investigating 
the photodynamics of Lil in the gas phase which is much 
less well characterized than that of Nal. The latter has 
been studied in detail, both experimentally using fem¬ 
tosecond spectroscopy, and theoretically using semiclassi- 
cal and quantum dynamical model calculations [8, 10, 22- 
27]. Preliminary results on Lil [9] suggest a similar pre¬ 
dissociation dynamics as for Nal. To clarify whether this 
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FIG. 1: Pump-probe scheme employed in the present experi¬ 
ments. The pump pulse excites a vibrational wave packet out 
of the groundstate X into the first excited state A. The probe 
pulse ionizes Lil by exciting into the groundstate of Lil + or 
into the Li + + I dissociative continuum near the outer turning 
point of the A-state. 

expectation is indeed correct, we perform one-color fem¬ 
tosecond pump-probe measurements in the UV spectral 
range. We corroborate the experimental results with 
high-level ab initio calculations to determine the poten¬ 
tial energy curves and non-adiabatic couplings from first 
principles. Based on this data, we model the photoioniza¬ 
tion dynamics of Lil, using both semiclassical arguments 
and full quantum dynamical simulations. 

Our paper is organized as follows. Section II in¬ 
troduces the electronic structure of Lil, outlining first 
the excitation scheme and discussing then the details 
of the electronic structure calculations. The experimen¬ 
tal setup and results are presented in Sec. Ill for gas 
phase measurements (Sec. IIIB) and droplet experiments 
(Sec. IIIC). We compare the experimental data to sim¬ 
ulation results in Sec. IV. Our findings are summarized 
in Sec. V. 


II. ELECTRONIC STRUCTURE OF Lil 
A. Excitation scheme 

The excitation and ionization scheme of Lil is depicted 
in Fig. 1 using the adiabatic potential energy curves 
obtained by ab initio calculations as described below. 
Analogously to Nal and other alkali halides, the struc¬ 
ture of adiabatic potential energy curves is determined 
by the low-lying Li + -I - ion-pair state which crosses the 
covalent states at intermediate internuclear distances. 
This causes the groundstate X to have a deep poten¬ 
tial well (Do = 28839cm -1 , R e = 2.39 A [11]) and the 


first excited state A to be predissociative. The cross¬ 
ing of the A-state potential with the ionic potential at 
Rx,a = 6.13 A [28] is avoided due to non-adiabatic as well 
as spin-orbit couplings with a coupling energy Vxa rang¬ 
ing between 686 cm -1 and 766 cm -1 according to the lit¬ 
erature [28-30]. A second curve crossing exists due to the 
intersection of the B -state potential with the ionic Li + - 
I - potential around R x ,b — 10 A where Vab ~ 5 cm -1 
according to our calculations. To the best of our knowl¬ 
edge, no literature data are available on the B-state. 

A vibrational wave packet, excited in the A-state by a 
first ultrashort laser pulse (left vertical arrow in Fig. 1), 
oscillates back and forth from the inner to the outer clas¬ 
sical turning points. Due to the avoided X-A curve cross¬ 
ing, part of the wave packet escapes outwards along the 
dissociative diabatic potential to form free Li + I frag¬ 
ments each time the wave packet comes near the crossing 
point. Since the ionization energy of atomic Li (5.39 eV) 
exceeds the photon energy (~ 4.6eV), only the bound 
excited molecule Lil* can be ionized by one-photon ab¬ 
sorption from a second probe pulse (right vertical arrow). 
The fraction of the wave packet which follows the adia¬ 
batic potential curve converging toward the ionic curve 
at large R can be ionized by photodetachment of the I - 
component, leaving behind Li + +I. Ionization is enhanced 
at the outer turning point of the wave packet oscillation 
around Ri = 12 A due to accumulation of wave packet 
amplitude at I?,;, which corresponds to a maximum clas¬ 
sical dwell time. At short laser wavelengths A < 270 nm 
the 13-state is also partly populated. This leads to fast 
dissociation along the diabatic 13-state potential curve. 


B. Ab initio calculations 

The electronic structure of the Lil molecule, briefly in¬ 
troduced in the previous section, is determined in state 
of the art ab initio calculations to allow for a quantitative 
comparison between theory and experiment. To ease cal¬ 
culations, the Lil molecule is treated as a system with ef¬ 
fectively 10 electrons. This is achieved by accounting for 
the core electrons of the iodide atom by a pseudopoten¬ 
tial. To account for relativistic effects, the ECP46MDF 
pseudopotential from the Stuttgart library [31] is em¬ 
ployed for the core electrons, leaving 7 valence electrons 
(5s 2 5p 5 ) in I which are explicitly treated in the calcu¬ 
lations. For the I valence electrons we adopt the stan¬ 
dard pseudopotential basis set [6s6p/4s4p] extended by 
[13d3f2g/5d3f2g\ functions taken from the quadrupole- 
quality basis set suggested in Ref. [32], For the Li atom 
the aug-cc-pVQZ basis set is employed [33]. The ab initio 
calculations of the Lil potentials described in the follow¬ 
ing are performed using the MOLPRO package [34]. 

In the first step we calculate the lowest nonrelativis- 
tic (S - A) states A 1 S+, (2) 1 S+, ^ 3 E+, and 3 II, 
see Fig. 2. These states asymptotically correlate with 
the groundstate atoms, Li( 2 S')+I( 2 P), and with the ionic 
asymptote, Li + ( 1 5 , )+I - ( 1 S'). The nonrelativistic singlet 
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FIG. 2: Calculated potential energy curves for the low-lying 
non-relativistic (S — A) states of the Lil molecule. The inset 
shows a close-up view of the potentials at short distance. 


(low-spin) states, X 1 E+, (2) 1 S + , and 1 II, have been 
calculated using the internally contracted multireference 
configuration interaction method with single and dou¬ 
ble excitations (MRCI) [35] with additional Davidson 
correction (MRCI+Q) [36]. The MRCI calculations are 
preceded by state-averaged complete-active-space calcu¬ 
lations (CASSCF) [37] to generate optimized orbitals and 
reference state energies for the MRCI runs. The multiref¬ 
erence method is essential here since we describe low-spin 
electronic states which asymptotically correlate with two 
open-shell atoms, 2 S + 2 P. The Davidson correction 
is highly important for the final results since it largely 
compensates for the lack of size-extensivity of the MRCI 
method. 

For the triplet states of Lil, 3 £ + and 3 II, and for the 
lowest states of the Lil + ion, 2 £ and 2 II, we could employ 
the open-shell spin-restricted coupled cluster method 
with single, double and non-iterative triple excitations, 
RCCSD(T) [38], which is a fully size-consistent method. 
The single reference coupled cluster method can be used 
here since these are the lowest high-spin states in a given 
symmetry. The long-range tail of the high-spin states, 
3 £+ and 3 n, can easily be obtained from the supermolec- 
ular RCCSD(T) calculations. However, the MRCI+Q 
results would not yield reliable interaction energies at 
long range and we need an alternative approach for low- 
spin states. For the ionic (2) 1 £ + state, the long-range 
part of the potential (R > 13 A) has been assumed to be 
pure — 1 /R Coulomb interaction. For the low-spin singlet 
states, X 1 £ + and 3 n, we have assumed the long-range 
interaction energy to be equal to the long-range part of 
3 £ + and 3 n states, respectively. This is a well justified 
approximation, since we know that van der Waals coef¬ 
ficients C n (n > 6) depend on the asymptotic atomic 
states but not on the resultant molecular spin; they thus 


are exactly the same for the low and high-spin states dis¬ 
sociating into the same atoms. The spin-orbit couplings 
between the nonrelativistic (S — A) states are calculated 
using the MRCI method with the spin-orbit Hamiltonian 
Hso in the Breit-Pauli approximation. The nonadiabatic 
radial coupling between the A 3 E and the (2) 3 E states is 
also calculated by employing the MRCI method. 

Finally we diagonalize the matrix H nre \ + Hso, where 
the elements of the diagonal H nre i matrix come from 
the MRCI+Q (low-spin states) or RCCSD(T)(high-spin 
states) calculations. After including spin-orbit couplings, 
the electronic states do not have any more a defined spin 
S or projection of the orbital angular momentum A on 
the molecular axis. The only conserved quantum number 
is now the projection of the total electronic angular mo¬ 
mentum (spin + orbital) on the molecular axis, denoted 
by fl. The H nre i+Hso matrix is block-diagonal, and each 
block can be labeled by a different quantum number Q. 
From the set of states we consider, A 1 E + , (2) 1 E + , 3 n, 
3 £ + , and 3 n, we obtain four blocks with fl = 0 + ,0 _ , 1 
and 2. Diagonalization of each of these blocks yields the 
adiabatic states with given O. For our purposes, the 
most interesting is the block with symmetry 0 + , which in¬ 
cludes the ground adiabatic state X0 + . Our model yields 
three 0 + states: A0 + , A0 + , and B 0 + , which asymptot¬ 
ically correlate to Li( 2 S’)+I( 2 P 3 / 2 ), Li( 2 S')+I( 2 Pi/ 2 ) and 
Li + ( 1 S’)+I~( 1 S'), respectively. 

The two 1 £ + states (or X0 + and A0 + after including 
the SO-coupling), strongly interact with each other near 
the avoided-crossing. To account for that, we calculate 
the nonadiabatic radial coupling between the two states, 
r(P), as a function of internuclear distance, P, 

r(P) = (A 1 E + |^ ? |(2) 1 E+). 

The coupling r(P) is then used to calculate the mixing 
angle y(P), 



In the next step, the mixing angle y(P) is employed to 
make a transformation from the adiabatic to the diabatic 
representation, 

V?(R) = F 2 ad (P)sin 2 7 (P)+F 1 ad (P)cos 2 7 (P), 
V d (R) = F ad (P)sin 2 7 (P)+F ad (P)cos 2 7 (P), 
F d (P) = (F 2 ad (P)-F 1 ad (P))sin7(P)cos7(P), 

where Vj ad (P) and V 2 d (R) are the adiabatic X 1 E + and 
(2) 3 E + potentials, respectively, F d (P) and V d (R) are 
the diabatic states, and Vj 2 (P) is the diabatic coupling 
potential. This diabatization procedure enables to elimi¬ 
nate the derivative with respect to P, r(P) • d/dR , from 
the equations which describe the nuclear motion [39]. In 
the dynamical calculations this derivative is replaced by 
the coupling potential V d 2 {R). It is highly advantageous 
since r(P) is a strongly varying function of P (approach¬ 
ing the Dirac delta form near the avoided-crossing), while 
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FIG. 3: (a) Dashed lines: Diabatic ab initio potential curves 
calculated in this work (see text for details); Solid lines: Adi¬ 
abatic potentials; Dotted lines: Potentials derived from spec¬ 
troscopy [28]. (b) Diabatic coupling functions. 


Vy 2 {R) is a smooth function of R. In Hund’s case (a) rep¬ 
resentation, the Hamiltonian governing the vibrational 
dynamics becomes 

(T + V^R) Vi d 2 (i?) <^ d |Hso| 3 n> \ 

# = V&{R) T + Vf (R) (F d |Hso| 3 n> 

V(F 1 d | J ffso| 3 n) (F 2 d |i?sol 3 n) t + iA + v 3 n (R)J 

(1) 

where T is the kinetic operator and A the energy splitting 
between the 2 P 3 / 2 and 2 Pi/ 2 levels of the iodide atom. 

Figure 3(a) illustrates the ab initio potentials (dashed 
lines) corresponding to the diagonal entries in Eq. (1). 
The coupling functions of the three coupled surfaces 
given by the off-diagonal elements in Eq. (1) are shown 
in Fig. 3(b). In case of the second excited state, the 
potential corresponds to |A + IA' n (P). Figure 3(a) ad¬ 
ditionally displays the adiabatic potentials obtained by 
diagonalizing Hamiltonian (1) as well as the potentials 
derived from high-resolution spectroscopy [28]. One no¬ 
tices a global very good agreement of the present theoret¬ 
ical results with the experimental curves. In particular, 
for the groundstate X0 + the present theoretical disso¬ 
ciation energy is 29543 cm -1 which compares well with 
28839 cm” 1 from Ref. [28]. The position of the minimum 
is 2.39 A from theory and 2.40 A from the experimental 
fit. For the H0 + state we can compare the T e parameter 
by determining the depth of the potential with respect to 
the ground vibrational level of the X-state. Our theoreti¬ 


cal value of T e = 30273 cm 1 for the A0 + state compares 
well with the experimental result T e = 30769 cm” 1 . 


III. EXPERIMENT 
A. Experimental setup 

The experiments are performed using a He nanodroplet 
beam apparatus which has been described in detail pre¬ 
viously [16, 40, 41]. The main modifications are a time- 
of-flight detector for measuring ion mass spectra as well 
as an amplified femtosecond laser system. The latter 
comprises a Ti:Sa-based oscillator (Tsunami by Spectra 
Physics) and a regenerative amplifier (Legend by Coher¬ 
ent) operated at 5 kHz repetition rate. Third harmonic 
radiation is generated by a home-built frequency conver¬ 
sion stage which consists of two BBO crystals, a calcite 
plate and a A/2-wave plate. The laser pulse characteris¬ 
tics are a pulse energy of up to 20 ^J, a pulse length of 
about 230 fs, and a wavelength tunability in the range 
A = 260-278 nrn. The laser pulses are split into identical 
pump- and probe pulses which are time-delayed using a 
mechanical delay line. 

An effusive beam of Lil is produced by heating a Lil 
sample to 380° C using a vapor cell which usually serves 
as a doping unit for the He droplet beam. In the experi¬ 
ments with Lil embedded in He droplets presented in the 
last section, the Lil cell is heated to 360° C for doping the 
droplets generated by a continuous expansion of He out 
of a 5 /iin nozzle by on average 0.8 Lil molecules. Photo¬ 
electron signals are recorded using a standard velocity- 
map imaging spectrometer setup [42, 43] in time-of-flight 
detection mode. 


B. Gas phase results 

Figure 4 displays the measured Li + (a) and Lil + (b) 
transient ion signals as a function of the delay time be¬ 
tween pump and probe laser pulses. The traces are ver¬ 
tically shifted from one another for the sake of visibility. 
The thin white lines depict the result of fitting the ex¬ 
perimental data to the function 

S(t) =S 0 + Sie” 4 /* 1 + S 2 e~ t/t2 sin (2 tt (t - t 0 )/T) . (2) 

It essentially consists of the sum of an exponentially 
damped oscillation (t 2 ) to model the wave packet dynam¬ 
ics in the bound state, and an exponentially decreasing 
offset (ti) to account for predissociation into unbound Li 
and I atoms which elude single-photon ionization. We in¬ 
terpret the presence of both Lil + and Li + product chan¬ 
nels by the population of both bound and unbound states 
close to the dissociation limit of Lil + , respectively, by 
photodetachment of I” in the Li + +I” ion pair. Exper¬ 
imentally, we find the ratio of Li + to Lil + ion yields to 
reach a maximum value of about 20 at the laser wave¬ 
length A = 270 nm. 
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Pump-probe delay [ps] 

FIG. 4: Measured pump-probe 7 Li + and LiI + -ion signal tran¬ 
sients recorded at various center wavelengths of the femtosec¬ 
ond laser. The smooth white lines depict fits to the data. 


The highest contrast of the wave packet oscillation in 
both Lil + and Li + transients is observed at the largest 
laser wavelength A ss 275 nrn. At that wavelength, 
damped oscillations up to about 8 ps are clearly visi¬ 
ble. A wave packet is excited deep in the potential well 
where dispersion due to the anharmonicity of the poten¬ 
tial curve is minimal. The subsequent reappearance of 
the oscillatory signal with maximum amplitude at about 
13 ps (“full revival”) results from partial rephasing of the 
wave packet and is determined by the anharmonicity of 
the potential [16, 44, 45]. Moreover, at that wavelength 
the excitation energy falls below the threshold for predis¬ 
sociation along the B-state. The question whether pre¬ 
dissociation via the B-state contributes to the reduction 
of signal lifetimes will be addressed below. 

For shorter wavelengths A < 275 nm, the oscillations 
quickly decay. The decreasing contrast is particularly 
pronounced for the Lil + signals, for which the periodic 
modulation turns into a fast irreversible drop of the signal 
within about 0.5 ps when exciting at A < 270 nm. This 
decay is assigned to direct dissociation via the repulsive 
diabatic potential of the B precursor state. The slow 
decline of the signal offset on the time scale of several 
ps barely changes as a function of A. In the following we 
concentrate on the discussion of the resulting fit values of 
the oscillation period, T, and the exponential decay time 
constant for the signal offset, t±. These quantities are 
depicted in Fig. 5 a) and b), respectively, as a function 
of A. 

As expected for a typical anharmonic molecular poten¬ 
tial, the oscillation periods are found to slightly decrease 
from about 0.85 to 0.6 ps for increasing A = 264-278 



Laser wave length X [nm] 

FIG. 5: Period (a) and exponential decay time constant (b) 
of the oscillating ion signals obtained from fits of the exper¬ 
imental data. The lines in (a) depict the classical oscillation 
period. The lines in (b) show the result of evaluating Landau- 
Zener’s formula for the diabatic transition probability at the 
avoided curve crossing between the X and the A-states. 

nm, see Fig. 5(a). Despite of the considerable scatter 
of the data points, the oscillation periods measured for 
the heavier 'Li + isotope systematically exceed those of 
6 Li + . No significant deviation of the oscillation periods 
measured for Li + and for Lil + are found within the ex¬ 
perimental error. 


C. Lil embedded in He nanodroplets 

In addition to measurements using effusive Lil, we 
report experimental results on Lil embedded in He 
droplets. The experimental setup is the same, except 
that the temperature of the vapor cell is slightly reduced 
and the beam of He nanodroplets is switched on. How¬ 
ever, the yield of Li + ions correlating with the He droplet 
beam as compared to the ion signals measured with ef¬ 
fusive Lil is reduced by about a factor 100. We attribute 
this dramatic reduction of the ion yield to the tendency 
of atomic ions to form tightly bound complexes (so called 
“snowballs”) inside the He droplets due to strong polar¬ 
ization forces [46]. 

Therefore, we switched to photoelectron detection 
which we find to be well-suited for recording the wave 
packet oscillations for free Lil molecules, as depicted in 
Fig. 6 (a). Indeed, the photoelectron signals for Lil em¬ 
bedded in He nanodroplets are found to reach about 30% 
of the signal level measured with effusive Lil (Fig. 6 (b)). 
However, no indications of a periodic modulation due to 
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Pump-probe delay [ps] 


FIG. 6: Pump-probe transients recorded by measuring the 
yield of 7 Li ions and of photoelectrons for Lil in an effusive 
beam (a) and for Lil embedded in He nanodroplets (b) at A = 
269,6 nm. The smooth lines in (a) are fits to the data. 


wave packet dynamics is observed. Note that we cannot 
exclude contributions to the electron signal from pho¬ 
toionization of contaminations in the He droplets, such 
as water picked up by the droplets from the residual gas. 
The lack of signal oscillations implies that the coherent 
vibrational wave packet motion is efficiently quenched by 
the interaction with the He environment. This contrasts 
previous experiments performed in our group where the 
droplets were doped with alkali metal dimers and trimers 
which reside in weakly bound dimple-like states at the 
surface of He droplets [13, 16, 47]. In those experiments, 
long-lived vibrational coherences were observed by de¬ 
tecting the neat dopant ions as a function of pump-probe 
delay. Thus, molecules immersed in the droplet interior 
appear to be subjected to much stronger perturbation of 
the vibrational dynamics than those sitting at the droplet 
surface. 


IV. MODELLING THE PREDISSOCIATION 
DYNAMICS 


at the temperature of the Lil sample (Tlu = 653 K), and 
H = 6.65 (5.74) amu is the reduced mass of 7 LiI ( 6 LiI). 
The resulting values (solid lines in Fig. 5 a)) reproduce 
the experimental ones reasonably well. This supports our 
interpretation of the wavelength dependence in Sec. IIIB 
in terms of the anharmonicity of the A-state potential. 

The decay constant of the signal offset, t \, is modeled 
using the standard semi-classical Landau-Zener formula, 


P(v x ) = exp 


ZkVxa 

HAVav x 


( 4 ) 


for the transition probability from one adiabatic poten¬ 
tial (A) to another ( X ) when a wave packet passes an 
avoided curve crossing. Here Vxa denotes the diabatic 
coupling energy between X and A-states, A Va is the 
slope of the difference potential at the curve crossing, 
and v x is the classical velocity at the crossing point in 
the center-of-mass system. Since in the present scheme 
the wave packet passes the X-A crossing point twice in 
each oscillation period, T, the total transition probabil¬ 
ity for non-adiabatic transition from the A to the A'-state 
per period amounts to 2P(1 — P ). This yields the time 
constant for the diabatic transition to the A-state, 


.LZ 
l i 


T 

2P(1 -P)' 


( 5 ) 


It is represented by lines in Fig. 5(b) for J = 30 and var¬ 
ious values of the coupling constant Vxa- The Landau- 
Zener model fits the experimental values well for Vxa = 
650(20) cm -1 , reproducing even the trend of a slightly 
increasing t-\ for larger values of A. 

To assess the contribution to the decay rate of the sec¬ 
ond A-B avoided curve crossing, we extend the Landau- 
Zener model to include the passage at both curve cross¬ 
ings. For the A-B avoided crossing we use R = 
10.6 A and Vab = 5 cm 1 as obtained for our ab ini¬ 
tio calculations. The resulting deviation of tf z with re¬ 
spect to the value obtained by accounting for only the 
X-A avoided crossing amounts to about 5%. Thus, the 
observed signal decay can be ascribed to predissociation 
predominantly via the lowest X-A avoided crossing. This 
is in agreement with expectations based on the nearly 
perfect diabatic behaviour of the reconstructed potential 
curves at the outer A-B crossing [28]. 


A. Semiclassical calculation 


B. Quantum dynamics calculations 


To confirm our assignment of the observed dynamics 
to wave packet propagation in the A-state of Lil we have 
simulated the vibrational motion classically for the A- 
state potential Va(R) augmented by the centrifugal term 

Vc{R) = h 2 J(J + l)/(2/rP 2 ). (3) 

Here, the rotational quantum number J = 30 is used to 
account for the thermal population of rotational states 


A more realistic modeling of the detected wave packet 
dynamics is achieved by treating the dynamics fully 
quantum mechanically including couplings between the 
relevant diabatic states and with the laser fields. In par¬ 
ticular our goal is to ascertain the role in the predisso¬ 
ciation dynamics of the second A-B curve crossing. To 
model the probe process in the presence of the laser field 
E(t), the Hamiltonian (1) is extended by adding the po¬ 
tential for the groundstate of the Lil + ion and the dipole 
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couplings d ■ E{t) between the two 1 S states of the Lil 
molecule and the ionic state. The dipole transition mo¬ 
ments d are assumed to be constant. 

The simulations employ the first three coupled diabatic 
Lil potentials and their coupling functions, as expressed 
by Hamiltonian (1), as well as the Lil + ionic groundstate 
potential. We assume that the observed wave packet dy¬ 
namics involves exclusively 0 + potentials. Initial pump¬ 
ing into the H = 1 states can be neglected due to the 
small transition dipole moment and the repulsive charac¬ 
ter of these electronic states. This approximation is jus¬ 
tified also by the fact that in spectroscopic observations 
of the transition X0 + —> A0 + there was no indication 
of perturbations coming from the f2 = 1 states [28]. To 
account for the thermally excited molecules, a rotational 
barrier parametrized by J is added to each potential ac¬ 
cording to Eq. (3). 

The pump process is simulated by defining a Gaussian 
shaped vibrational wave packet in the first excited state 
at t = 0. The peak position and width are determined 
by mapping the groundstate vibrational wave function in 
the X-state onto the A-state potential according to the 
spectral intensity distribution of the laser pulses used in 
the experiment and the Va{R ) — Vx(R) difference po¬ 
tential. For simplicity, only the vibrational groundstate 
is assumed to be populated prior to the excitation step. 
Thus, the calculations include propagation of the initial 
wave packet and interaction with the time delayed probe 
pulse. Consequently, the population transferred to the 
ionic state via single-photon absorption is calculated as 
a function of delay time. Since the rotating wave approx¬ 
imation cannot be applied in these calculations, the time 
grid is required to have a resolution of 97 as to accurately 
sample the oscillating electric field of the ultraviolet laser 
pulses. The emission of photoelectrons upon excitation 
into the ionic state is not taken into account. 

Due to the predissociative character of the first two ex¬ 
cited states, artificial reflections of the wave packet at the 
edge of the numerical grid occur which are suppressed by 
adding complex absorbing potentials. Following Ref. [48], 
we use a quadratic complex potential, 

V (r) = —Z 77 • 8(R - R c ) ■ (R - R c ) 2 , ( 6 ) 

for the dissociative first two excited states, for which re¬ 
flection plays a role. Equation ( 6 ) employs the Heaviside 
step function 6 and is characterized by the parameters 77 
(strength of the potential) and R c (onset of the poten¬ 
tial) . 77 and R c need to be chosen such that all dissociat¬ 
ing fractions of the wave function at a sufficiently large 
distance are absorbed. 

The time-dependent Schrodinger equation is solved by 
expanding the time evolution operator in Newton poly¬ 
nomials [49]. The pulse parameters are chosen to match 
the experiment in that the pulse length (full width at 
half maximum, FWHM) t p = 230 fs is determined by 
an autocorrelation measurement whereas the intensity of 
the pulses in the focal region is estimated to / « 2 x 10 11 
W/cm 2 based on the measured beam diameter. For a 



FIG. 7: Experimental and simulated pump-probe transient 
7 Li+ ion yield at the center wavelength A = 266.7 nm. The 
smooth lines are fits of expression ( 2 ). 

detailed comparison with the experiment, the resulting 
transient signals for rotational states up to J = 70 have 
to be thermally averaged according to the Boltzmann 
distribution of Lil at T = 673 K, 

{S)=Y j Sj9jPj/Y,9jPj- ( 7 ) 

j j 

Here, Pj = exp (—Ej/ksT), Ej = hcB 0 J( J + 1) denotes 
the Boltzmann factor, gj = 2J + 1 is the degeneracy of 
the J-th rotational level, and Sj is the simulated delay 
time dependent population in the ionic state for given 
J. A comparison of the pump-probe traces at the wave¬ 
length A = 266.7 nm of the experiment and a simulation 
for which every second rotational state between 6 < J < 
70 is taken into account for simplification is shown in 
Fig. 7. The traces show a good agreement in terms of 
the period and damping of the oscillations as well as the 
exponential decay. Also, the vibrational revival occurs at 
a similar time of roughly 10 ps. 

For further simplification, the calculations for the 
wavelengths investigated in the experiment are per¬ 
formed with a single state, J = 30, only. The simu¬ 
lated pump-probe transients are analyzed in the same 
way as the experimental data by fitting model Eq. (2). 
The resulting simulated oscillation periods T and decay 
constants t\ are shown in Fig. 8 where they are compared 
to the experimental values. The agreement between the 
simulations and the experiment is fairly good, confirming 
the accuracy of the ab initio potentials and couplings. 
The simulated values of T shows the same dependence 
on the wavelength. The decay constant t\ becomes min¬ 
imal for A ~ 269 nm as in the experiment. However, 
slight differences between simulation and experiment are 
present. In particular T lies systematically below the 
measured values. This can be attributed to neglecting 
any initial vibrational excitation of Lil in the quantum 
dynamics simulations. The results of simulations, where 
















FIG. 8: Fitted experimental (filled symbols) and simulated 
(open symbols) oscillation periods (a) and exponential decay 
time constants (b). Solid lines represent the classical oscilla¬ 
tion periods (a) and the Landau-Zener model (b) as in Fig. 5. 

full thermal averaging over rotational states is performed, 
are depicted in Fig. 8 as well. Thermal averaging does 
not significantly improve the agreement with the experi¬ 
ment which justifies the simplified approach using a single 
state, J = 30. 


V. SUMMARY 

We have studied the predissociation dynamics of Lil 
in a one-color femtosecond pump-probe experiment for 
molecules in the gas phase and embedded in He nan¬ 
odroplets. The wavelength was varied in the range 
A = 260 -278 nm. We have corroborated the experimental 
observations by full quantum dynamics simulations based 
on newly calculated ab initio potential energy curves and 
non-adiabatic couplings. Vibrational wave packet oscilla¬ 
tions in the predissociative A-state are observed in both 
ion and electron yields for gas phase Lil molecules. The 
exponential decay of the signal with increasing pump- 
probe delay times is interpreted in terms of predissocia¬ 
tion due to diabatic coupling between the excited A-state 
and the ground X-state. The good agreement of the ex¬ 
perimental and simulated pump-probe transients confirm 
the accuracy of the calculated potentials and coupling 
energies. Predissociation is found to occur nearly exclu¬ 
sively via the X-A avoided curve crossing. Since it is a 
single crossing that determines the predissociation, the 
decay constant as a function of laser wavelength can also 
be reproduced by a semiclassical Landau-Zener model for 


the diabatic transition. The reason that the dynamics is 
dominated by the X-A avoided crossing is due to the 
strong diabatic coupling between these states, found to 
amount to 650(20) cm -1 . This is much larger than the 
coupling between the A-state and the H-state of only 
5(2) cm" 1 . 

Analogous measurements with Lil embedded in He 
nanodroplets result in photoion signals that are strongly 
suppressed compared to the gas phase measurements. 
This is attributed to complex formation inside the 
droplets upon ionization. While the photoelectron sig¬ 
nals have comparable amplitudes to those measured for 
free Lil, no wave packet oscillations are observed. This 
may be explained by a strong coupling of the excited 
Lil to the He droplet which induces fast damping of the 
coherent vibrational wave packet motion. 

Our observation is in contrast to femtosecond pump- 
probe experiments with diatomic alkali metal molecules 
attached to He nanodroplets for which the vibrational 
motion is well resolved [13-17]. The difference is ra¬ 
tionalized by the different coupling of the molecules to 
the droplet environment. It results from different loca¬ 
tions, on the droplet surface for the weakly coupled alkali 
dimers whereas molecules immersed in the droplet inte¬ 
rior such as Lil are subjected to strong coupling. For 
molecules inside the droplets, a similarly strong effect of 
the He environment on the molecular dynamics as ob¬ 
served here has also been reported for recent impulsive 
alignment experiments [21], The droplets were found to 
alter the rotational period and suppress alignment re¬ 
vivals. 

These findings together call for a more comprehensive 
study of the ultrafast dynamics of molecules immersed 
in a He quantum fluid. For our pump-probe experiment, 
refined detection schemes such as ion imaging techniques 
would allow to unequivocally confirm suppression of co¬ 
herent wave packet motion. More generally, it would 
also be interesting to see whether the regime of interme¬ 
diate coupling between molecule and He droplet can be 
explored by using different molecular species, e. g. alkali- 
earth alkaline mixtures [50]. This would provide much 
needed input for the theoretical modeling of the inter¬ 
action of molecules with He droplets and thus be crucial 
for developing a rigorous understanding of quantum fluid 
environments. 
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